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CHAPTER I 
IKTRODUCTIOK 



1.1 General Statement of the Problem 

When a nuclear weapon explodes, a high pressure wave develops 
and moves avay^from the point of the e:q>losion, traveling over the ground 
svurface and through the ground. A shallow- bur led arch subjected to this 
pressure resx>onds in varioxxs ways, depending on the orientation of the 
structure relative to the direction of travel of the blast wave, the parame- 
ters of the blast loading, and the properties of the soil suirounding the 
arch. If the arch is oriented with its longitudinal axis i>erpendicxilar to 
the direction of travel of the blast wave, the windward side of the arch 
is subjected to the pressure a few milliseconds before the leeward side. 

This unsymmetrical sequence of loading produces flexural stresses euid de- 
flection in the arch. The deflection of the arch into the svirrounding soil 
gives rise to an unsymmetrical 'passive earth pressure loading. 

An understanding of how these loadings occur and their effects on 
the arch is necessajry for the proper consideration of the design of the 
arch rib. 

1.2 Purpose of this Thesis 

When the arch rib is subjected to the above described loads it 
responds in a configuration which is a function of an infinite number of 
shapes or modes of vibration, since the arch rib has distributed mass and 
elasticity. An infinitely large number of coordinates would be required 
to specify the position of the deflected arch. The participation of each 
mode depends on the spatial distribution of the loads along the arch rib 
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as well as the time variation of the dynamic load, in relation to the natural 
frequencies of vibration of each mode. 

The distortion of the arch rib is a function of the stiffness of 
the rib, the paasive earth pressvtre loading resulting from the stiffness of 
the soil surrounding the arch, the distribution of the soil and dynamic 
loads, the ti^ variation of the dynamic load, and the support condition of 
the abutments. It is convenient to consider the earth pressxare loading due 
to the stiffness of the soil as an integral part of the resistance of the 
arch. For flexible arches, the stiffness of the soil is greatly in excess 
of the stiffness of the arch rib, consequently the properties of the arch 
rib may be neglected. Henceforth, unless specifically designated otherwise, 
the resistemce of the arch rib means the resistance given to the arch by 
the soil. 

Design of structures for dynamic loading is usually based on the 
jarinciple of conservation of energy. When the dynamic load acts on the 
structural element, the element deflects. Work is done on the structure by 
the load at any point sis that point deflects. Due to the speed of the load- 
ing, kinetic energy of motion is given to the mass of the element, and, as 
the element deflects, strain energy is stored. The strain energy may involve 
both elastic and plastic action. A general statement of the principle of 
, conservation of energy is : 

work done = kinetic energy + strain energy 

To evaluate this equation the displacement of any point on the 
sttmcture \inder the load must be known. The strain energy is the sum of 
ell the strain energy in t he element, and may be caused by conpressive, 
flexural, torsional and shear stresses in the element. The kinetic energy 
is the energy of translation or rotation of all the masses in the element. 
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The evaluation of these quantities would be exceedingly difficult 
if no more than the arch rib were involved. Evaluation of the energy 
stored in the soil is a hopeless task, impossible for practical purposes. 
Fortunately, the principal modes of response of the arch rib are a sym- 
metrical mode corresponding to the concession stresses and an anti symmetrical 
mode corresponding to the flexural stresses. This allovs a considerable 
sln^liflcation to be made in the treatment of the problem. 

A single -degree-of- freedom system consists of a lun5>e<i mass 
restrained by a weightless spring. When this system is subjected to a load, 
only one coordinate is required to fully define the position of the mass. 

By assuming that the anti symmetrical mode of response of the arch rib is 
defined by the radial deflection of a point at the quarterpoint of the arch, 
^diich is the point of maximum deflection of the eirch, the entire deflected 
shape of the arch may be defined with one coordinate, and the arch may be 
considered a single-degree-of-freedom system. In essence this neglects any 
tangential movements of the arch rib associated with the radial deflection. 

After representing the action of the arch rib as that of a single- 
degree-of-freedom system, dynamic equivalence relations should be established 
between the actxrnl and idealized masses, spring stiffnesses, and loads acting 
on the two systems. However, for the purposes of this study, this refinement 
is omitted. Instead, the following £issunq>tions are made: 

(1) The mass of the single-degree-of-freedom system is equal to 
the unit mass of the arch, which is obtained by uniformly distributing along 
the arch rib the weight of the soil which creates the passive earth pressvire. 

( 2 ) The resistance function of the single-degree-of-freedom 
system is equal to the unit resistance given the arch by the earth pressure. 




-■m 



k 



( 5 ) The forcing function on the single-degree-of -freedom system 
is equal to the unit pressure of the forcing function acting on the arch rib. 

Once the mass, resistance function, and forcing function acting 
on the single-degree-of- freedom system have been defined, the response of 
the system may be determined by the well-knovn differential equation of 
motion of such system: 

na + kx = p(t) 



where Sd* 
kx 
m 

P(t) 



the second derivative of the displacement with respect 
to time 

the resistance of the system as a function of the 
displacement 

mass of the system 

the externally applied force as a fxuiction of time 



The main purpose of this thesis is to investigate the response of 
various arch ribs under a limited range of variables when acted upon by blast 
pressure. In order to do this, it is necessary to relate the parameters of 
the single-degree-of-fTeedom system to those of the actual structure under 
investigation. 



1.5 Method of Approach 

The blast pressure acting on the surface of the groxmd causes 
vertical and horizontal pressures in the soil. These pressures are a 
function of the maximum surface pressure, and are time-dependent upon a 
pressure wave which passes throiogh the soil. Pressure- time relations of 
these pressures are formulated. 

These pressures are assumed to be principal stresses, the resultant 
of \hich acts radially on the arch. The tangential component is omitted. 

The total radial load idiich results from these pressures is designated as 
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a load profile Pj.(<P>t), which is then represented as a distributed load by 



the series 



I 



sin 



nwcp 



where sin Rd <p 



$ - central angle of the arch 
R » arch radius 



An analogoxis, uniformly-distributed, antisyrametrical, time- 
dependent load, Pg^(t), is then determined. This load has the same general 



effect on the entire haunch of the arch as the load 



I 



sin 



nK 9 



The pressvire-tiitte diagram for p (t) and the numerical relation between p (t) 

8> SL 

and the dynamic load ) P„(t) sin — are established by equating their 






effects at variotis times as though they both acted statically. The ansilogous 
antisymmetrical load p„(t) defines the forcing function acting on the arch 
and also on the single-degree-of -freedom system. 

The unit, horizontal, passive earth pressure acting on the arch 
is determined by Rankine *s theory for passive earth pressure. The unit 
vertical pressure acting on the arch is foxind by statics. These two pressures 
are assumed to be principal stresses, and their radial resultant is assumed 
to resist the movement of the arch. The tangentieil consonant is omitted. 

The total radial soil load is designated as a load profile p^(9), and is 



represented as a distributed load by the series 

e 



Z p sin 



n« <p 



where 






2 

e R 



J p,(i>) 



sin Rd cp 

U 



An analogous, uniformly-distributed, antisymmetrical load which 
has the same general effect on the entire haunch of the arch as the load 
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Z nJt cp 

sin - is determined by statics. The passive earth pressure is 
assumed to reach its maxi m\m value •tdien a point located at Q/k on the arch 
rib deflects radistlly an amovint equal to O.OOIR. In addition^ the resistance 
diagram of this pressxire is assvimed to be elasto-plastic . The analogous 
anti symmetrical load defines the vmit resistance of the arch rib and also 
the resistance function for the single-degree-of-freedom system. 

The veight of the soil that causes the passive earth pressure is 
determined by Rankine's theory and by statics. This weight is assumed to 
be uniformly distributed along the arch rib and defines the unit mass of the 
arch rib and also the mass of the single-degree-of-freedom system. 

? The response of the equivalent single-degree-of-freedom system is 

determined by a method developed by Dr. N. M. Newmark for blast resistant 
design^ Reference ( 8) . 
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CHAPTER II 

general discussion OT ^ VARIABLES 

The variables of the problem may be divided into three general 
groups corresponding to the blast loading, the earth pressure loading, and 
the configuration of the arch rib. 

2.1 Blast Pressure Loading 

The blast forces acting on the arch rib are conceived to act as a 
pressure wave passing through the soil. Figure 2.1. The direction and 
velocity of the propagation of this wave in the soil is ass\mied to be related 
to the shock velocity of the surface pressure wave, U, and to the seismic 
velocity of the soil sxuTovmding the arch, C . The magnitude of C is 
assumed to be less than U. 

Within the stressed-soll zone, the vertical pressure p^(t) is 
assumed to be equal to the surface overpressure P^(t), but to have a dif- 
ferent pressure-time curve to allow for the retardation effects on the shock 
transmitted through the soil. The horizontal pressure Pjj(t) is assumed to 
be proportional to the vertical pressure. 

The time of arrival t^ of the pressure wave at any point in the 

soil, or on the arch, is determined by the propagation velocity of the 

wave. The time of rise t of the vertical pressure to its maximum value P 

r m 

is assumed to be a function of the vertical transit time t^ of the pressure 
wave. The rate of decay of the earth pressure wave is aissumed to be the 
same as that for the surface pressxnre. These relations are shown 
qualitatively in Figxare 2.2. 

Study of Figure 2.1 indicates that as the shock velocity of the 
blast wave increeises relative to the seismic velocity of the soil, the 
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unsyrametrical loading decreaises, since the pressure wave approaches the 
arch DK>re nearly vertically. For the usual range of these variables the 
earth pressure wave qualitatively approaches the arch as shown in Figure 2.1. 

It sho\ild be noted that this pressure wave is considered to have 
been induced, solely by the surface pressxire. The direct ground shock wave 
associated withT a nuclear weapon e^qploslon is not considered in this thesis. 
Reference (l) states that for pressures of less than 300 psi and for depths 
of less than fifty feet the direct ground shock may be ignored. 

The arch is assumed to be acted upon by the radial resultant of 
the vertical and hori^ntal pressvires. If the ratio of the horizontal 
pressure to the vertical pressure is as great as 1:1^ the radial pressiare 
is equal to the vertical pressure. Moreover, the total radial force in- 
creases as the ratio of the horizontal to the vertical pressure increases. 

The effect of the rise time may be seen in Figure 2.2. The 
pressure wave arrives at some point on the arch, Figure 2.1, and the 
pressure increases to its maximum value in some time t^. At some adjacent 
point <pg the pressure is at a different level, depending on the relative 
arriveLl and rise times. Reference (2) recommends that the rise time be 
taken equal to one-half the vertical transit time. 

The following range of variables will be studied in this paper: 

= 100 and 200 psi with U = 5CKX) and 4000 fps 

Cg ■ 1000 and 2000 fps 

t ■ ■fc* 1/2 t. 

IT V V 

Pjj/Py ® 0*25> 0.50, 0.75 and 1.00 

2.2. Passive Earth Pressure Loading 

The unsymmetrlcal blast load on the arch causes the arch to deflect 
into the soil on the leeward side. This deflection brings into action the 



9 



passive earth pressure of the soil as veil as the inertia of the mass of 
soil associated with this pressure, since the blast load acts dynamically. 

This pressure contributes materially to the resistance of the arch and must 
be considered. 

While precise evalxiation of the passive earth pressure under static 
conditions is C9H^)lex, evaluation under dynamic conditions is equally, or 
even more, con5>lex. The passive earth pressure can vary in the length of 
the arch due either to the properties of the soil or to the construction 
technique used in placing the backfill. Rankine’s theory has been chosen as 
the sinq>lest method for evaluating this peissive earth pressure. 

In the application of Rankine's formulas. Reference ( 5 )^ the unit 
weight of the soil is assumed to be 120 pounds per cubic foot and the angle 
of internal friction is assumed to be thirty degrees. 

2.3 Structures Considered 

The blast loading on arches with central euagles of I80 and 120 
degrees and radii of 150 and ^00 inches is investigated to determine the 
effect of the soil variables, arch geometry, and relative rise time on the 
arch loading, and to determine a conventionalized dynamic loading function. 
Subsequently investigated is the response of arches with central angles 
of 180 and 120 degrees and radii of 200 to 600 inches to this conventionalized 
dynamic loading function. 



5 
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Time 



Figure 2.2 Vertical Pressure-Time Diagram for 9 
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CHAPTER III 



GENERAL RELATIONS COMMON TO EACH LOADING FUNCTION 



To predict the effect of the loads on the arch rib, the loads must 

be defined in a convenient manner. 

Certain general relations may be derived for defining either the 

blast load or the earth pressure load. These general relations concern the 

transformation of stress vith direction, the method used to determine the 

value of the coefficients p (t) and p for the distributed radial loeid 

•*^n' ' ^n 

functions ^ — and ^ p^ sin which represent the actual 

load on the arch, and the basis for assuming analogous loads to represent 
these distributed loads. 

For purposes of dynamic analysis, the arch is Idealized to a 

Z i \ nJt Q) 

p^(t) sin — g — ■ 

is replaced by a unifonnly-distributed, antisymmetrical, time -dependent load 
p^(t), and the distributed earth pressure load ^ p^ sin is replaced 

by a uniformly-distributed antisymmetrical load p . These antisymmetrical 

6l 

loads act outwardly on one half of the arch, and they act inwardly on the 
other half. 

The direction of the analogous earth load is opposite to that of 
the analogous blast load. These antisymmetrical loads define the effects 
of the actual loads on the arch haunch. 

The values of these analogous loads are found by the application 
of the principles of statics. The value of the analogous load for the soil 
pressure may be found by equating the moment produced by the distributed load 



I 



p^ sin — at a point 0/4 to those of the suialogous antisymmetrical losid 



•n 



p^ at the same point. The replacement of the dynamic loaui^p^(t) sin ^ 
with an a nal ogous load requires that a basic assiunption be made. It is 
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assumed that if tvo generally similar static loads produce equal static 
effects on a structure at a specified point, then the same tvo loads applied 
dynamically with similar pressure-time diagrams of equal duration vould 
cause equal dynamic effects at the same point. Of course the dynamic 
effects of the loads are not defined. The point considered is 0/4. 

To delJermine the analogous loads, the moments caused by the anti- 

symmetrlcal loads Pg^(t) and p^ are equated vith the mcaaents caused by the 

Z , . nit <P n>t cp 

p^(t) sin - and \ p^^ sin ^ respectively. In 

addition, the effect of a uniform compressive load on the central third of 

the arch coupled with a uniform outward load on the outer thirds of the arch 

is determined. This load provides an analogous load for study of the crown 

of the arch. 



5*1 Stress Transformation with Direction 

Let it be assumed that a horizontal unit stress Pj^ and a vertical 
vinit stress p^ are acting on a segment Rdcp of a semicircular arch rib located 
at the angle <p measured from the right abutment. Assuming that p^ and Pj^ 
are principal stresses, the radial pressure Pj.(<P) at <p may be found by the 
well-known method for determining the relations between stresses at a point 
on different planes passing throxigh the point. Reference ( 4) . 

Rd9 may be eissximed to be a plane sxarface for an infinitesimal 
segment. Hie vertical projection of the plane Rd9 is Rd<p cos <P; the hori- 
zontal projection is Rd<p sin 9 . Equating the sum of the radial forces to 
zero yields the following: 

O p 

Py(9) Rd9 = p^ Rd9 cos 9 + p^ Rd9 sin 9 

Using the trigonometric identities 

2 

cos 9 = ^1 + cos 29 ) 

p 

sin 9 * i(l - cos 29 ) 
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axid cancelling Rdcp from both sides of Equation (a), 

Pj.(*P) “ ^1 + cos 2qp) Pjj + 5(1 - cos 2<p) 

= + Py) + KPh - Py) cos 2(f> (1) 

The ratios of Pj^/P^ ^ ^ studied are 0.25^ O. 5 O, 0.75 and 1*0. 

In addition, these pressures are time -dependent, so Pj.(<p) must be a time- 
dependent pressure and is therefore expressed as p^(9>t) when the pressures 
considered arise from the blast pressure anting on the surface of the ground. 



Let Pjj(t) * Cp^(t). Then Equation (l) becomes: 

Pj.(9,t) = C p^(t) + p^(t) (1 - C) sin^cp (2) 

This equation is evaluated for various values of C as follows: 

For C » 0.25, pJ<P,t) = ^.25 + 0.75 sin^cp) P^(t) (5) 

For C = 0.50, p^(9,t) =(0.50 + 0.50 sin^9) P^(t) (4) 

For C - 0.75, Pj.( 9 ,t) =(0.75 + 0.25 sin^cp) P^(t) ( 5 ) 

For C = 1.00, p^(9>t) = p^(t) (6) 



Equations (5) through (6) express the time -dependent unit radial 
pressure acting on the arch rib at any point 9 for the various assumed 
ratios of the horizontal pressure P^('t) 'to the vertical pressure p^(t). 
Evaluation of these equations is set forth in Table I. 

5.2 Distributed Load F unctions 

Any radial load profile p^(9), such as that shown in Figure 5*1^ 



when Eicting upon a segment of a circular arch rib with a central angle 0, 
may be represented as a distributed load by a trigonometric series in the 
following manner: 
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Pj.(9) = sin ^ + Pg sin 



2 n 9 



n«9 



^ 4- — - + Sin 



(7) 



The coefficient p^ is obtained by multiplying both sides of 



/ \ nJi 9 

Equation (7) by sin — 7 — and integrating over the length of the arch rib. 



J Pj.(«P) sin ^ m q> = J p^ sin ^ si 



n«9 

sin — ^ Rd9 

a 



I 



+ ^ ?2 sin ^ sin ^ Rd9 

o 



+ + 






. 2 n «9 

sin Rd 9 + — 



(8) 



L 



e 



xjjjcT) nJtO) 

The value of the integral j sin -- ^ - c in Rd 9 is zero provided 
that m and n are integers and m does not equal n. Equation ( 8 ) then yields: 



2n>»9 



sin H2 Rdcp . R [ I - ] 

6 



i p„ R0 



’ k f '’r('f) ^ 



(9) 



If P^(<P) can be ejqaressed analytically^ the value of p^ may be 
found by integration. If p^(9) is an irregular load profile for ^diich 
there is no convenient analytical expression, the integral may be evaluated 



by numerical methods since the integred 



/Pr(9) 



nX9 

sin Rd 9 represents the 

U 



area luider the load profile multiplied by sin 



n?tp 



The sign convention adopted throughout this thesis is that a 
positive load acts inwardly and a negative load acts outwardly. The 
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n 

harmonic components of loading p sin shovn in Figiare 5*1 are in 

n u 

qualitative agreement vith the original load p^(<P)j therefore the sign of 
the coefficient Pg is definitely negative. The sign of the coefficient p^ 
has been taken as negative, although no positive statement can be made that 
this is the case for the original load p^(<P) in Figure 5*1* 

Included in Figxare 5*1 are tvo cases of uniform reidial loading. 

p sin on the crovn and haunch 

n e 

of the arch may be represented by such uniformly distributed loads. 

If the load profile p^(?) is time-dependent, it is ejqjressed as 
p^(<p,t). The procedure described above may be used to determine the co- 
efficient p^(t), with the qualification that P^Ct) must be evaluated at a 
specified time. 

Reference (5) contains a full development of the methods of 
harmonic analysis as applied to structural engineering problems. 



5.5* Basis for the Dstermination of Analogous loads 

As stated at the beginning of this chapter, the analogous loads 
p (t) and p are derived by equating their moments to the sum of the moments 

€L £L 

of the distributed loads ^ and > respectively. 

The time dependency of p (t) is not considered in this section since the 
load ^ pjt) sin is assumed to be acting statically at the instant 

the analogous load p (t) is ev£d.uated. 

Moment equations for all load profiles in Figure 5*1 are found 
in this section after general equations for the right reaction of any two- 
hinged semicircular arch are established. Since the blast loeid and soil 
load are e^qpressed in terms of a general function ^ p^ 
of moment equations suffices for both analogous loads. 



sin — ~ , one set 

U 
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5.5.1. General Equations for the Right Reactions of C ircular-Segment Arches 

Since the circular- segment two-hinged arch. Figure 5*2, is one 
degree indeterminate, one equation based on Castigliano's theorem. 

Reference (U), is used for the determination of the reactions. The horizontal 
tlrrust at the right reaction is chosen as the redundant. It is assumed that 
the arch rib may* be treated as a thin ring, making it unnecessary to consider 
the theory of curved beams. The effects of rib shortening also are not 
considered. 

The following geometrical relations are derived in accordance with 
Figure 5*2: 



Figinre 5* 2a 



Figure 5* 2b 



Figure 5* 2c 



Figure 5* 2d 



BC = 2R sin^ 

AB = BC cos|-{ 0 - cp) = 2R sin^ coc^ 0-9) 
AC = BC sin^ 0 - 9) = 2R sin^ sin^ 0-9) 
BD = 2R sin^0 

DF = 2R sin^e - ^) 

EF = 2R sin^0 - ^) cos^0 - 0) 

= moment arm of P about D 

FG = R sin (9 - p) 

= moment arm of P about G 

AB = moment arm of V about C 
AC = moment arm of H about C 
DB s moment arm of V about 0 



In Figure 5. 2d, the total strain energy U in the arch rib due to 
the moment M caused by the radial load P acting at a point on the rib 
located by p may be e:q)ressed as 
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Since the abutments are assumed not to move, the contribution of the right 
reaction H to the strain energy is zero and may be expressed eis 










0 



(a) 



Two equations are required to express the moment in the arch rib; 
one equation for angle cp less than or equal to p, and one for cp equal to or 
greater than p. For 0 < <? 5 p, 



= V(AB) - H(AC) 

= 2VR sin ^ cos^ 0 - *?) 

3m, 

= -2R sin^ sin^0 - cp) 



- 2HR sin^ sin^e - cp) 



(b) 

(c) 



For P < cp < e. 



Mg = V(AB) - H(AC) - P(FG) (d) 

= 2VR sin^ cos^e - cp) - 2HR sin^ sin|{ 0-cp) - PR sin(«P-0) 

^Mg 






= -2R sin^ sin^P - q>) 



(e) 



Substituting Equations (b) through (e) into Equation (a), 



^ Rdcp+^Mg ^ Rd«p 



^Mg 



= r [ 2VR sin^ cos^ 0-cp) - 2HR sin^ sin^ 0-cp) ] 
'^o 

[ -2R sin^ sin^0-9) ] Rd<p + 

O0 

/ [ 2VR sin^ cos^ 0-cp) - 2HR sin^ sin^ 0-cp) 

'^p 

- PR sin(cp-p) ][-2R sin^ sin^0-cp)] Rd9 



(f) 



Solving Equation (f) for H yields 
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H = 



4 J VR^ sin^ I sin(-~) cos(-^)d9 

r® 

-2 / PR^ sln(cp-p) sin | sin(-^)d<p 
4 r 5 sin2 I sin2(^) dq> 



(g) 



Referring to Figure J.2, the sum of the moments equated about 
the left reeiction gives 

P(EF) 2PR sin(-^) cos(-^) p 



DB 



2R sii^ 



2 sin I 



( 10 ) 



Substituting Equation (lO) into Equation (g), 



J sin( 
4 si 



i-n(e-fi) 



. e 

sin^ 



i - sine 



- i cosej + l^cos I - ^ cos(| - p) 



H = 



- I cos(|^ - P) -1- (^) sin(| * P) j ]■ 



e ^ « 2/ ex . e e 

^ + e cos (^) - 5 sin ^ cos ^ 



( 11 ) 



Equations (lO) and (ll) define the right reaction of any circular- segment 
two-hinged arch acted upon by a radial load P located at any angle p. 



5.5*2. Moment Equations for I80 Degree Arches 

For a semicircular arch, the central angle e is n , and 
Equation (3JL) reduces to 

H = ^ (sin p + (p - «) cos p] (12) 

Using Equation (12), the horizontal right reaction may now be 

n 

found for the harmonic load con^jonent sin . 



ir^i 

i „ 








19 



JRp 



For P = p, sin Rdcp = sincp Rdcp^ since 0 = «, 
i. C7 X 



, r « PiR 

H = i / Pj^ sin cplsincp + (cp-«) coscpl Rdcp = 



For P = -Pg sin Rdcp , 



H = -Pg sin2cp[sincp + (cp-«) coscp] Rd9 = ^ PgR 



For P = -Pj sin jcp Rdcp^ 

if* X 

H = — J -Pj sin5cp[sincp + (<P-k) cos<p] Rdcp » g p^R 

n JCp 

The loading configurations p sin — ~ are shovn in Figure 5.1. 

n C7 

It is also necessary to investigate the effect of the uniform 
antisyrmnetrical load p^ and the uniform symmetrical load p^ shown in 
Figure 5»1. Fc5t the uniform antisymmetrical load, 
a ^ 

H = J ^ -p [sincp + (cp-n) coscp] Rd<p + ^ p[sincp + (cp-a)coscp)Rdcp = pR 



For the uniform symmetrical loeui. 



a 2a 

H = ^-p[sincp + (cp-a) coscp] Rdcp + p[sincp + (cp-a)coscp] 



Rdcp 



^ ^a P Islo'P t (^-*) coscp] Rdcp = j -n/ 5 pR 



Equation (lO) is used to find the right vertical reaction. 
With 0 = a. Equation (lO) reduces to V = 2 R sincp. 
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P ^ p 

For P = sin ^ Rdq), V = ^ J sin^cp Rdcp = ^j^«R 



For P = -Pg sin ^ Rd<p, V = - Pg sin2<p sin<p Rdqp 



= 0 



r ** 

For P = -p^ sin ^ Rdq>, V = - 5 / p, sinjcp sincp Rdq) 



= 0 



For the uniform anti symmetrical load, 

1C 

n 









-p sincp Rd<p + i L p sincp Rdcp = 0 



For the uniform syrametriceil load 

« 2H 

V = ^ -p sin Rd + z P sinCp Rdcp 



- r 



p sincp Rdcp = 0 



The value of zero for the laat four loads vas to be expected from 
the symmetry of the arch and the loading configurations. It is interesting 
to note that the two-hinged semicircular arch is statically determinate for 
the last four loadings. 

All the above values for the right reactions for a l 8 o° arch are 
tabulated in Table II. 

The moment at any point on the arch in Figure 5*1 located by <P 
due to a distributed load p, sin ^ Rd 0 acting at p is 

i. C7 

dM = -p^ sin 3 [R sin (cp- 0 ) ] Rd^ 



The total moment at cp is 




r 

I 
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r9 

M = -H(AC) + V(AB) - J ?! sin(9-p) 1 Rd^ 

» - 2R sin I sin(^^j + 2R sin | cos(^) j 

+ 6in<p sin^ ^ - cosqp ( | - ) j 

^ X - (fl-2<p) cos9 - 5 sin<pj 



u Pi 



A similar procedure is used to find the moment at <p due to the 
other loading configurations in Figure J.l. ®ie results are tabxalated in 
Table III. 

5 . 5 . 5 . Moment Equations for 120 Degree Arches 

When the central angle of the arch is 120 degrees. Equation (U) 

reduces to 

sin(~ - 



) 5 " fi a , 2« 1 2x1 

. 5 31 . 3. - i cos 3-j 






+ COS I - I cos(| - p) - ^ cos (a-p) + sin(| - p)| 

® ” a 2a 2 a I a a ^ 

s- + -s- cos T " 5 sin -r cos •=■ 

3 5 5 5 5 

P j^(cos0 + ^ 6inp)( I - ^ «) + i -(^ + cosp 

+ (^ - sinp + ^ p cos P - P sinp J 

1 ^ 

2 



= P 



1.85979 " 2.70578 cosp + 0.16665 sin0 
+ !• 59550 p cos p - 0.91990 P sinpj 



(15) 
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Equation (ij) may be used to find the right horizontal reaction 
of an arch with a 120 degree central angle acted upon by a concentrated 

j(Cp 

radial load P located at any angle p. The reaction for sin Rdcp is 
fo\ind as follows; 



H = 




Pi sin 



9 J^l. 



85979 - 2.70578 coscp + 0.16665 Sincp 



+ !• 59330 9cos<p - 0.91990 <P8in<p 



Rd<p = 0. 53076 p^R 



The horizontal reactions for the remaining harmonic loads and the 

uniform loads are found in the same manner, a n d are set forth in Table II. 

Equation (lO) is used to find the vertical reaction. 

«cp 

For Pj^ sin ~ Rd9 , 

2n 

V = J ^ \_ — ~ — « j ^1 2 ^1^ 

sin ^ 

The remaining vertical reactions have been fo\md, and are stated 
in Table II. 

For a segment Rdp at p loaded with p^ sin Rd P , the moment 
at 9 in Figure 3*1 is 



dM = -p^ sin ^ P[R sin(9-P) 1 Rdp 
and the total moment at 9 is 

M = -0.55076 p^rUc] + 0.6p^ RtAB] 



- J^'Pi si® I PtRsin{(p-p)) Hap 



r — - <P 1 

= -0.53076 p^R 2R sin I sin(-2 — ^ ) 



] 



2n 

T- - 9 



+ 0.6 Pj^R j^2R sin ^ cos(-^ ^ 

“ Pi sin 2 P iR sin(<p-p) ] Rdp 

= (0.78499 - 1.35965 sincp - 0.78499 cos9 + 0.8 sin p^R' 
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The moment equations for the remaining load configurations shovm 
in Figure 5 .I for an arch with a 120 degree central angle have been foxxnd 
and are contained in Table III. 



5 . 3 Procedure in Determining Analogous Loads 

The moment equations tabxilated in Table III have been evaluated 

for every ten degrees. The results are tabulated in Table IV and shown 

in Figures 3*3 3 **^« These values are for static loads. They do not 

indicate the dynamic moment which may exist at any time in the arch. 

Inspection of Figures 3*3 3 *^ shows that the effect of the 

analogous antisymmetrical load is almost identical to that of the load 
2«Q> 

Pg sin — ^ . If Pg has a value of one, then the value of the antisymmetrical 
load must be slightly less than one if the effect of each is to be the same. 

However, the antisymmetrical load must also include the effect of the loads 

jcp 3»p ^ njcp 

p^ sin ~ and p^ sin . Therefore, the moment that \ p^ sin causes 

at a point Q/h is first determined, and subsequently a value of the anti- 

symmetrical loaul is found viiich will cause the same moment at 0/4. The 

point 0/4 is selected since the msymmetrical blast load produces the maximum 

effect at this point. The time must remain constant for each n\imerical 

evaluation of p (t), as described below. 

The actual blast load on the arch is p^(9, t). At any given time 

the load profile p^(9>t) can be determined and analyzed to give a distributed 

load ^ p^(t) sin . By ansuming that ^ acts statically, 

its effects may be found at 0/4. An analogous load can be found •vdiich 

causes the same effect. By repeating this operation at varioxis times, the 

pressure-time diagram of p (t) can be determined. This load must then be 

Si 



€q>pUed dynamically to determine its actual effect on the arch. 
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In determining the static load of the earth, the analogous load 
is found by equating the moments at 0/4 of the distributed static load 

sin — ^ with the moment of the anti symmetrical load. Since the earth 
pressure is not a dynamic load, only one evaluation is needed. 

Further inspection of Figures 5*3 and J.4 shows that the effect of 

Z nJtcp 

p^ sin — ~ on the central third of the arch may be 
represented by the symmetrical load which acts in compression on the central 
third of the arch and in tension on the outer thirds. This load is 



designated p 



If the distributed load is a time-dependent load, the same 



procedure is used in finding the press\are-time diagram for p (t) as that 

8 

used for p ( t) . 

£l 
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TABLE I 

COEFFICIENTS FOR RESOLVING VERTICAL AND HORIZONTAL PRESSURE INTO 

RADIAL PRESSURE 



Angle 


Radial 


pressure Pj,(9>t) in terms of vertical 


pressure 






Assumed 


ration p^/p^ 






0.2^ 


0.50 


0.75 


1.00 


0 


0.250 


0.500 


0.750 


1.00 


10 


0.273 


0.515 


0.758 


1.00 


20 


0.338 


0.559 


0.779 


1.00 


30 


0.458 


0.625 


0.815 


1.00 


ko 


0.560 


0.707 


0.855 


1.00 


50 


0.690 


0.794 


0.897 


1.00 


6o 


0.815 


0.875 


0.938 


1.00 


70 


0.915 


0.942 


0.971 


1.00 


8o 


0.978 


0.985 


0.993 


1.00 


90 


1.000 


1.000 


1.000 


1.00 


100 


0.978 


0.985 


0.995 


1.00 


no 


0.913 


0.942 


0.971 


1.00 


120 


0.815 


0.875 


0.958 


1.00 


130 


0.690 


0.794 


0.897 


1.00 


lltO 


0.560 


0.707 


0.853 


1.00 


150 


0.458 


0.625 


0.815 


1.00 


l6o 


0.558 


0.559 


0.779 


1.00 


170 


0.273 


0.515 


0.758 


1.00 


180 


0.250 


0.500 


0.750 


1.00 



REACTIONS OF CIRCUIAR-SBGMENT ARCHES WITH H AND V IN TERMS OF pR 
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(a) Radial Load Profile p^(9) and Harmonic Load Comgponents 
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(b) Uniform Symmetrical and Anti symmetrical Loadings 



Figure 5 .I Various Load Configurations 



50 






DE H To FO 
E F iso Tangent Line 




Figure 5*2 Geometry of a Circular Arch Segment 
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Figure 5*5 Moment vs Angle 9 for l80 Degree Centreil Angle 
Two-Hinged Arch 
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Figure 5.4 Moment vs Angle cp for 120 Degree Central Angle Tvo-Hinged Arch 
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CHAPTER IV 



BLAST PRESSURE LQADIHG 

4.1 Pressure-Time Relations at a Point 

As stated in the general discussion of the variables. Section 2.1, 
the blast pressure acting on the svurface of the ground above the point of 



interest causes a vertical pressure at that point. Accon5>anying this 
vertical pressure is a horizontal pressvu:e of varying magnitude, depending 
on the soil. The horizontal pressure may be as lov as fifteen percent in 
loose, dry soil, fifty percent in moist soil, and one hvindred percent in 
saturated soil. Reference ( l) . One of the purposes of this thesis is to 
stvidy^the effects of this variation. 



equal to the seismic velocity of the soil. At the point of interest, the 
time for the pressure to increase to a maximum is a function of the vertical 
transit line. The effect of varying the rise time will be studied. 



may be e:q>ressed in terms of the vertical pressure. Vertical presstire-time 
relations at a point on the arch in terms of surface pressure are there- 
fore required to define the total radial load profile p^(9, t) actliag on 
the arch. The equations below are derived from the geometry of Figure 4.1. 



The vertical pressxire is propagated downward with a velocity 



It has been shown in Section J.l that the horizontal press\n:e 



The velocity of the wave propagation is; 




( 14 ) 



where 



V 



velocity of propagation of the earth pressure wave 



Cg = seismic velocity of the soil 
U = surface bleist wave velocity 



34 



The time required for the pressure wave to arrive at any point 9 
on the arch rib is: 



where 



is: 



t = R - R 6in(cp - p) 
a V 

9 = angle locating a point on the arch rib 
R = radius of the arch 

t = time meaisured after initiaO. contact of the pressure 
wave with the arch rib 

C 

p = arc tan -rr 



( 15 ) 



The time retjuired for the pressure at 9 to rise to its peak 



R - R sin9 + H 

c 



R - R sin9 + H 

t = ^ S 

r 2Cg 



for t = t . 
r u 



for t = 5 ta 



(16) 



^ere 






= depth of cover on the crown 
t^ = rise time of the pressure 

t = vertical transit time of the pressvure wave to 9 



t 



t 



R - R sin9 + H 
c 

C 



s 

The vertical pressure-time variation at a point on the arch rib 
is shown qualitatively in Figure 2.2. The pressxare-time diagram is assumed 
to be a straight line in the rise part of the diagram and to have an ex- 
ponential decay rate after reaching its peak vaJLue The total inq?ulse 

of the force pulse at any point must be the same as the total inq)\ilse of 
the svirface force pulse. This means that as the rise time increases^ the 
decay rate must also Increase in order to maintain the same impvilse in 
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the force piilse. It is assumed in this study that the decay rate at the 
shallow points in the soil considered in this thesis is the same as the 
decay rate of the surface force pulse. This is a reasonable ass\ui 5 )tion 
since the maximum unsymmetrical loading occurs within the first few milli- 
seconds after the blast wave begins to act on the arch. 

The de&ay rate for a side-on overpressure having a positive 
phase duration t^ may be obtained flrom Reference (6). The decay rates for 
the two overpressure levels stxidied in this thesis are shown in Figure k»2. 
For the 100 psi overpressure, t^ is assiuned to be I 5 OO milliseconds, \diile 
for 200 psi overpresstire, t^ is assumed to be l61»0 milliseconds. The 
weapon yield is 1 ME. 

Hie complete vertical pressure-time diagram for any point can be 
determined only by numerical evaluation of the pressure on that point at 
various times. A sample calculation of a pressure-time diagram follows. 

Assume that a I 80 degree central angle arch with a radius of 
twelve and one-half feet is subjected to a 200 psi oveipressure blast wave 
originating from a 1 MP weapon. The duration of the positive phase of the 
loading is l 6 ^ milliseconds. The arch is buried in soil having a seismic 
velocity of 1000 fps, with zero depth of cover on the crown. The time of 
rise of the earth pressure at a point in the soil is assumed to equal the 
vertical transit time of the earth pressure wave. The ratio of the 
horizontal to vertical pressure is assumed to be 0 . 25 . 

Reference (7) gives the following equation for the shock velocity 
U of the surface blast wave: 

U = cjl + 6p/7P^)2 

where c^ = ambient sound velocity ahead of the shock 
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p = peak overpressure behind the shock 

F s asiblent pressure ahead of the shock 

X 

.*. U = 1150 (1 + 1200/105)2 . 4025 fFs 

uc 

By Equation (l4), v = — — - = ■ . . = 972 fps 

y Cg + y 1000^ X 4025^ 

By Equation (15), B » arctan (1000/4025) = 15® 58' 

t a [1 - 6in(cp - 15°58')] , seconds 

= 12.86 [1 - sin(cp - 15®58'3 , milliseconds 

*r “ ^5o ^ ^ * seconds 

= 12.5 (1 - sincp), milliseconds 

40*^. Then 

t = 12.86 [1 - sin(40® - 15°58’)] = 7.2 ms 

Gl 

t^ = 12.5 (1 - sin 40®) = 4.5 ms 

Therefore at cp = 4o degrees, the earth pressvire wave arrives 
7.2 ms after the wave initially contacts the arch rib at <p = 90® + 15® 58' . 
The vertical pressxire rises at a rate of 200/4.5 = ^5 psi per ms. At 8 ms, 
the pressxore is (8 - 7*2)45 = 56 psi. At 10 ms, the pressure is 
(10 - 7.2)45 = 126 psi. At 11.7 ms, the pressure reaches a maximum value 
ecjial to the surface overpressoare, or 200 psi. Thereafter the pressure 
decays at a rate assumed to be equal to the decay rate for the surface 
overpressure as shown in Figure 4.2. Inspection of Figure 4.2 shows that 
the decay may be represented ewjciarately by a straight line down to a value 
of A P/^P|jj = 0.65 and ^t/t^ = ~ 0.02. The pressure decrement is then 



By Equation (16) 



Let 9 = 
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0.55(200) = 70 pel, and the time Inteirval is 0.02(l64o) = ~ 55®s* So the 

pressure decays at a rate of 70/55 = 2 psi i>er ms. Therefore at 12 ms, the 

pressure is 200 - 2(l2 - 11.7) = 199 psi* At llf ms, the pressxrre is 

(199 - 4) = 195 psi. This procedure is continued until a time greater than 

(55 + 11) ms is reached. Then the time interval is expressed as a fraction 

of t, and AP/AP 'is read from Figure 4.2. The conrolete pressxjre-time 
d m 

relation for <P = ^ degrees is recorded in Appendix B, Table B5. 

To obtain the con5>lete vertical pressure-time relations for the 
entire arch, the above procedure is acccMZ5>lished for evejry ten degrees on 
the arch. Table B5, Appendix B, shovs the verticeG. pressure-time variation 
for the entire arch up to 200 milliseconds. The tmsymmetrical loading is 
readily seen. For exaisple, at ten milliseconds the vindvard side of the 
arch has pressvire acting on its entirety, •vhereas the leevard side is not 
fully loaded until after sixteen milliseconds. The unsymmetrical loading 
shifts from the vindvard side to the leevard side during the time the blast 
pressure acts. This effect may be seen at thirty milliseconds. The 
pressure on the vindvard side, cp greater than ninety degrees, averages some 
eight potinds less than the press\ire on the leevard side. The effect of 
the decay rate may be seen by con^>aring the pressxnre at the crovn vith 
that at the abutments. For example, at tventy-four milliseconds the pressxire 
at the crovn has reached its maximum and decayed to approximately seventy- 
five percent of the maximum surface pressure, vhile the presstire at the 
abutment is approximately equal to the maximum surface pressvjre. 

4.2. Harmonic Analysis of RadieGl Blast Load vith Time a Constant 

Based on the vertical pressures tabulated in Table B5# a radial 
load profile Pj,(9^t) can be determined for any instant of time given in 



» 
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that table. After p (9/t) has been defined, the distributed load 

r 



I 



/ \ n*P 

Pjj(t) sin — ~ can be evaluated. A sample calculation, for t = 10 milli- 



seconds, is given in Tables V and VI. 

In Table V, the data in column "p^( t) " is obtained fipom Table B5, 
•vdiile that in colimm ”p^ coefficient" is derived from Table I for Pjj/p^ “ 
0.25. The data in column "p^(<P>t)" of Table V is derived by multiplying 
Py(b) by the p^ coefficient. The resulting radial pressure load profile 
Pj.(9/‘t) is shovn graphically in Figure The coefficients P^(t) of 

Equation (9) are determined by nvunerical integration. This is done by 
multiplying the load ordinate at each ten degree interval by ^ sin — ~ , 



adding the resxilts and multiplying the sum by 



8R 



The resulting dis- 



tributed radial loading^ ' with t = 10 milliseconds, is 

shovn as a dotted line in Figure 4.5 • 

Table VI contains the harmonic analysis of the same blast vave 
acting on a 150 inch radius arch with a central angle of 120 degrees. The 
radial pressure p^(9/t) in Table V from 50 degrees through 150 degrees is 
used in Table VI from 0 degrees through 120 degrees. The resxilting dis- 






tributed radial loading \ p (t) sin is shovn in Figure 4.5. 






^n 



0 



Inspection of the data in colxann "Pj.(9, t)". Table VI, reveals that 
there is a radial pressure of ten or more pounds acting at all points on 
the arch rib. Therefore, a uniform all-around pressure of ten pounds may 
be subtracted from the radial press\ire before the harmonic analysis is made, 
since it is desired to resolve harmonically only the radial pressxare causing 
flexure in the arch rib. The column "p^ adjusted" contadns the remaining 
radial pressure. 

The coefficients Pjj(t) of Equation (9) are found by numerical 
integration. This is done by multiplying the load ordinate at each ten 

degree intei*val by ™ sin 2^ , add i ng the results, and multiplying by . 

^ OR 
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The resulting distributed load, * vith t = 10 milll- 

seconds, is shovn as a dotted line in Figure k.^. 



1 



It is assumed that the first three terms of the series 

Q JCCp 

sin — adequately represent the loading on the arch for the purposes 



of this thesis. 

The radial load on the l80 degree arch at a time t = 10 milli- 
seconds is: 

p^(<P,t) = l66 sin9 - 15 sin2<p - 26 sinjcp 

Since the time is constant, the load may be considered to be a static load. 
An analogous losid may be determined 'vAiich will have the same effect on a 

given point of the arch. Let the point be located at 9 = ^5 degrees. 

Z D 

p^ sin all cavise 

negative moment at this point. Using the coefficients in Table IV, the 
moment of the load is 



M= [l66( -0.02262) + 13( -0.55555) + 26( -0.08857) 1 
= [-5*75 - 4.55 - 2.50 ] = -10.58r^ 



A viniformly distributed anti symmetrical load acting in on the 
windward half of the arch and out on the leeward half produces a moment of 

P 

0.4l422p R at 9 = 45 degrees. Therefore 

a. 



10.58r^ 

0:4i422R^ 



25 psi for a time t « 10 ms 



However, the actual load on the arch is time-dependent and consequently p 

£l 

must be a time -dependent load, p (t). To define the pressure- time diagram 

£l 

for p (t), the load p (9^t) on the arch is harmonically analyzed at variotis 

£L X* 

times. 



The analogous load for the crown may be found in a similar manner. 



The moment caused by the above load at the crown is 



M = [166(0.055140) + 26(0.125) ] 

2 

=» 9*15 R inch pounds 

A symmetrical load p eicting in compression on the central third 

of the arch and in tension on the outer thirds of the arch produces a 

2 

aoment of 0.15l»'7 P R at the crovn 

8 

2 

• P«("t) = -5 = 59 psi> t = 10 ms 

® 0.15l^7R 

4.5 Blast Loading Approximations 

Appendix B contains the results of numerous harmonic analyses of 
the blast load on an arch arising from a 1 ME nuclear veapon. The analyses 
vere made and analogous loads foimd as described in the preceding section. 
These analyses vere made in order to determine the effects of the variables 
discussed in Chapter II on the loading the arch receives, and to provide a 
basis for idealizing Ihe blast load on the structxire vithout the necessity 
for actually harmonically analyzing the blast load. 

The blast load must be defined in relation to the effect being 
studied. In this thesis, the effect studied is the deflection of the arch 
due to flexxiral stresses, the area of primary interest being the haunch of 
the arch. The effect of the unsymmetrical blast load on the haunch of the 
arch is most conveniently related to an analogous uniformly-distributed 
load \diich acts invardly on one half of the arch and outwardly on the other 
half. The effect of the blast load on the crown is conveniently related 
to a uniform load acting in ccmipression on Idie central third of the arch 
and in tension on the outer lliirds. Values of such loads for various 



assumed conditions are shown in Tables BU and B12 
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After the blast load has been related to the specific effect 
being stiidied, a load-time cvarve must be established to shov the variation 
of the load vLth time. It viDJL be noted in Tables Bll and B12 that the 
analogous loads are for a specified time only, and that they vary appreciably 
with time. Appendix B contains a few typical analogous load- time cui*ves. 
Similar plots were made of all analogous loads. 

Study of these load-time curves revealed that the peak values of 
the analogous loads are some fairly uniform function of the maximum surfew^e 
overpressure for all assumed values of the variables. In addition, the 
shapes of these load-time cinrves are approximately initially peaked tri- 
angular load pulses, \diich decay rapidly to a relatively small value and 
then continue to decay at a much slower rate. The duration of the rapid 
decay phase of the loading corresponds approximately to the time required 
for the pressure to engulf the arch. The duration of the slow decay phase 
of the loetding is relatively much longer. 

The peak values of the analogous loads found in Appendix B are 
shown in Table VII and plotted in Figures 4.4 and 4 .^. The data in 
Appendix B euid Figures 4.4 and 4.5 may be qualitatively summarized as 
follows; 

(a) For ai5y given overpresstire, the peak value of the 

analogous load for the haunch increases as the value of 
the seismic velocity of the soil increases. This is ex- 
plained by the fact that a higher seismic velocity allows 
the pressure to penetrate the soil to a greater extent 
and thus a larger segment of the arch on the windward side 



is loaded before the leeward side is loaded 



For any given overpressiare, the peak vgilue of the 
analogovis load for the crovn decreases as the seismic 
velocity of the soil increases. !IMs is eacplained by 
considering the effect of the load profile p (9#t) on 

pcp ^ 

the integral / ^ P„(9^t) sin Rd<p . When C has 

sj r C7 s 

a relatively high value in con 5 >arison with U, a greater 

portion of the segment from to is located on the 

windward side of the arch and the values of the sine 

function of the angle 9 are both positive and negative. 

When C has a lower value in comparison with U, the 
s 

load profile p^(9/t) is located more nearly in the central 
third of the arch, 'rfaere the value of the sine function 
of the angle 9 is negative. This Increases the numeri- 
cal value of the analogous crown load. 

The higher ratios of Pj^/P^ increase the peak value 
of the analogous load for the haunch. This is explained 
by the fact that the area of the total radial load profile 
p^(9>t) is greater for the higher ratios, and thus the 
total unsymmetrical load is greater. 

The Pjj/Py ratio has little effect on the peak value 
of the analogous crown load. The peak value of the crown 
load is primarily a function of the location of the load, 
i.e., the central third of the arch. In this area, the 
value of the coefficient for determining the resultant 
of the horizontal and vertical forces is approriraately 
one for all values of Pj^/p^. Consequently, changing the 
ratio has little effect on the radiel load. 



(e) The assumed rise time has relatively little effect 
on the peak value cf the analogous loads. The peak 
values of both loads are reached very early in the load- 
ing of the arch vhen the pressure vave has not penetrated 
the soil to any extent. For shallow depths, the rise 
time is' practically equal tofee rise time of the stirface 
pressure, tdiich is assumed to be instantaneous. 

(f) Although the higher ratios of Pj^/Py cause a higher 
peak value of the analogous haunch load, they actually 
cause a less severe dynamic load effect. If the ratio 
is high, the pressiore, after the pressure wave has en- 
gulfed the £LTch, is more nearly a vuaiform compression, 
causing less flexural action in the arch. On the other 
hand, if the ratio is low, the tmlformity of the compres- 
sion is reduced, increasing the flexureO. action in the 
arch. For the higher ratios of Pj^/Py it is found that 
the analogous loads change direction, thereby creating 

a negative phase of loading. The peak value of the 
pressure in this negative phase is small, but could pos- 
sibly cause an outward movement of the crown of the arch 
since it has a relatively long duration. 

U^.3.1 l8o Degree Arches 

Since the blast loading may act over all, or over part, of the 
arch segment, one approximation to the peak analogovis loads may be obtained 
by considering the blast loading acting on the arch in such a way that the 
maximum value is given to ecich conq)onent of loading in the series 
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Z p (t) sin . The ratio of the horizontal pressure to the vertical 

pressure which causes the nost severe dynamic loading is 0.25^ therefore 
this value is used in the following derivations. The effect of the rise 
time of the pressure and the time of arrival of the pressvire wave is not 
considered. The vertical pressiore at all points on the arch is assvuned to 
he equal to ihe vmit radial pressure at any point is, hy Equation (5), 

Pj.(9,t) = (0.25 + 0.75 sin^cp) 

Using Equation (9), 

Pn(^) ^ J ^ 

2P 2 

maximum Pj^(t) = I (0.25 + 0*75 sin'^cp) sin9 Rdcp = 0.955 Pjjj 

'Jq “ 

maxinaan Pg(t) = 1 ^^ (0.25 + 0.75 sin‘^q>) sin2Cp Rd<p = -O. 39 Q 

2 

2P 

maximum p^(t) = (o*25 + 0.75 sin'^cp) sin39 rd9 = -O.4o8 P^ 

5 

The maxi nn an effect of the blast loading which may occxar, vinder an assxaned 
ratio of Pj^/P^ = 0.25, can never exceed the effect of the sum of maximum 
values of the harmonic components fovuid in this manner. The value of the 
individual maximum values found, and the effect of the sxan, are always 
reduced when the seismic velocity of the soil, the shock velocity of the 
surface blast wave, and the pressure-time relationships are considered. 

If these con^wnents acted statically and separately, they would 
produce the following moments at the leeward haimch emd crown: 



M = - j^0.025(0.955Pnj)R^ +0.5^3(0.598Pm)R^ + 0.088(0.iK)aP^)R^ 

= -0.190 P 
m 

= 0.055(0.955P^)R^ + 0.125(o.4oaPj^)R^ = 0.085 P^R^ 



The peak value of^ the analogous haunch load is 



Pa(t) 



0.190 P R^ 
^ m 

0.4llt R^ 



0.46 P 

m 



The peak value of the analogoxis crovm load is 

0.084 P 

" - 0.15 1 1 7- “ 

These values of "the analogous loads can never occur, and, for 
want of a better name, they are designated as the maximax analogous loads. 
They are shown in Figure 4.4 so "that they may be con5)ared with the analogous 
loads foijnd by actxial harmonic analysis. Such coiqiarisons indicates that 
these approximations are conservative. The following reductions are 
arbitrarily made: 



Maximum p (t) = 0.37 P 
■^a' ' m 


(17) 


Maximum p (t) = 0.45 P 

S 211 


(18) 



Equations (17) and (18) are plotted in Figtire 4.4. 

To obtain an idea of the magnitude of these analogous loads, they 
may be related to a uniformly distributed load acting on a sin^ily supported 
beam. Inspection of Figure 5*5 shows that the moment caused by the analogous 
load p on the haunch of the arch is very similar to the single beam 
moment. The moment at 45 degrees caused by the a n alogous load is 

M = 0.57(0.41422) PjjR^ 



U6 



Assume that the simple beam ha« a length equal to the length of half the 
arch^ or L = -^ • The moment is vL^/8. Equating this to the moment in 
the arch, it is found that: 



V = 



8 X 907 X 0.4ll4-55 P R' 



= 0.5 P 
m 



A similar analogy may be drawn for the analogous load acting on 

the crown, with the exception that the length of the simple beam is now 

SR 

one-third the length of the arch, or . The moment in the arch at 

90 degrees due to the analogous load is 






M = 0.45(0.1547) PjjR' 



8m 



8 X 0.45 X (0.1547) PjjR" 



(f) 



“ 0-5 

The above maxLmax loads occur within approximately two to four 
milliseconds after the loading begins to act on the arch. This fact is 
evident by inspection of the data in Appendix B. These loads decrease 
rapidly since p^(t) becomes practically zero \dien the pressure wave has 
engulfed the arch, and p^(t) also decreases. However, Pj^(t) and p^(t) are 
still acting. It is now eissumed that the pressure has con5>letely engulfed 
the arch. Since the pressure has engulfed the arch, the all-arovuid viniform 
pressure may be subtracted from the e:qiresslon for p^.(9>t). Then 

Pj.( 9 ,t) = 0.75 sin^<p Pj^ 



f 




t 




This radial pressure may be resolved into harmonic con^wnents. It is 
found that 



Pl(t) 

Pj(t) 



2P 



m 

« R 



2 P 




J (0.75 sinS) 

(0.75 sin^q)) 




sinqj Rdcp = O .657 P 

in 

sin39 Rd<p = -0.127 P 

m 



The laoment in the haunch at *^5 degrees caused by these two 
comopoi^Qbs is 

= - j^0.025(0.637PjR^ + 0.088(0. 127Pjj^)R^j = -0.026 P^^R^ 

The haunch analogous load is 



pjt) - 



0.026P 

m 

"■ o.ii.iir' 



0.06 p 

m 



( 19 ) 



The coniponents p^(t) and p^(t) act for the total duration of the 
surface blast pressure. Reference ( 2 ) states that the duration of the 
positive phase of the side-on overpressure may be considered as 



vhere 



W 



O.to W (20) 

Pm 

veai>on size in MP 

the effective duration in seconds of a tria ngu l ar pulse 
having the same impulse as the actual shock. 



The total duration of p^^(t) and p^(t) is considered to be t^. 

The time required for Pg(t) to become practically zero is found 
by vising Equations (15) and (16), with the angle 9 * 0 . 
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This time is defined as: 

tj *» t + t = duration of initial triangular inipulse 
d a r 



R + R sin p j. R 



V 






•with t = 



it. 



The time t^ defines t for "the load P^vt) = 0.06 



4.5.2 120 Degree Arches 

Approximations of the blast load ewting on a 120 degree arch may 
be derived in the same manner as that used in finding the approximations 
for "the 180 degree arch, as described in Section 4.3.1. It is necessary 'to 
modify Equation (5) by substituting (cp + 50) for since in Equation (5) 

9 is meas\ired ftrom a horizon-tal radial line Joining the arch abu'tment with 
the cen'ter of the circular segment. Equation (5) becomes 



Pj.(<P,t) = [0.25 + 0.75 8 in^(<p + 30)1 

= 9 + ^ sin9 cos9 + ^ cos^9 j P 

Let p^(t) = Pjj.(9,t) sin ^ Rd9 



Then the maxlnium possible values of PQ('t) are 



^ 2 « 

Pp — / 

Maximum p, (t) = ^ ( r ^ sin^cp + sin9 cos? 

i-*R ^ ® 






= 1.103 p 

m 



+ ^ cos^9) sin ^ 9 Rd9 



2P ~ ^ 

Maximum p«(t) = ^ sin^qp + sin9 coscp 

£«r Jn 



/3 



r«R Jt 



+ ^ cos^cp ) sin5<? Rdqp 



= -0.505 P 



m 

4 k 



Maximum p,(t) = ^-2 / ^ ( E ^ sin^9 + sinqp cos9 

^ ~KR '^2Jt ^ ^ 

+ ^ cos^qp ) sin I 9 Rdqp 



_-KR ''2K 

^ 9 



= -0.417 P 



m 
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The harmonic con^ponents Pjj(4) sin — ^ vould produce moment at the 
haunch and crovn as shovn belov^ if the blast load is considered to act 
statically. The moment coefficients are obtained from Table IV, for cp = 50 
and 60 degrees. 



«h= - 



0.009(l.l05Pj^)R^ + 0.125(0.505P^)R^ + 0.050(0. 417P^)R^ 



J 



0.094 P R" 

m 



M^ = 0.0l4(l.l05P^j)R^ + 0.051(0.417 Pj^)R^ 



= 0.057 PjR^ 

Using the moment coefficient for 9 = 50 and 60 degrees in 
Table IV, the analogous loads are found to be 

0.094 P R^ 

Maximum p (t) = - 

0.154TR^ 

» 0.61 P 

m 



50 



Maximum p (t) = 
8 



0.057 P 

m 



0.06556r^ 



= 0.56 P 



m 



These values are plotted in Figure 4.5 so that they may be 



compared vith the'actual analogous loads found by harmonic analysis^ as 
given in Appendix B. Since they are conservative, the following arbitrary 
reductions are made; 



These equations are plotted in Figure 4. 5* 

These loads may be related to a sinple beam load of 0.5 P^^^ by the 
same method as that described in section 4.5.1 for a I 80 degree arch. 

After the pressure wave has engulfed the arch, there is a uniform 
all-aroxnad pressure acting which must be subtracted from the radial 
pressure in order to determine the loading which causes flexural action. 
This uniform pressure is 



The remainiiag radieO. pressure vAiich is to be resolved harmonically is 



Maximum Pj^(t) = 0.45 P^^^ 

Maximum p (t) = 0.4l P 
■*^s m 



( 21 ) 



( 22 ) 





Then 
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2P 

Pl(l) = J' ^ ^ Bin^9 + 8in<p co8<P + ^ cos^tp) 8in|<P Rd<p 



0.51*6 P 



PjCt) 



m 



2P 



+ ^ sin^9 + ^ sin<p co8<P + ^ cos^<p)8in^ Rdq) 



= -0.020 P 



m 



These loads vould produce the follovlng moment at the haunch and 
crovn of the arch if they acted statically: 

= - |^0.009(0.546p^)R^ + 0.050(0. 020PjR^j 



= -0.006 P r' 
m 



= 0.0l4(0.5l*6Pjjj)R^ + 0.051(0.020 Pjjj)R^ 



= 0.009 P R 

HI 



The analogous einti symmetrical load for the haunch is 

0.006 P R^ 



Pjt) = 



m 



O.I55R 



0.04 p 



m 



The analogous load for the crown is 



(25) 



O.OO9P R 

P.(t) 



o . o 6556 r ^ 



0.14 P 



m 



(24) 
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The time t that defines the analogous loads in Equations (25) 
and (24) is the time required for the pressure vave to engulf the arch a n d 
rise to its maxiinum value at the abutments. This time is expressed by 
Equations (15) and (16). The angle <P in both of these equations is thirty 
degrees . 



4.5.5 Summary 

The blast load acting on the arch may be defined by nxuaerically 
evalxiating at specific times the coefficient "the distributed load 

p^(t) sin load- time curve for the analogous loads is obtained 



z 






by equating the static effects of the distributed load ^ ^ 

those of the suialogous lo€ids. This is the procedure employed in Appendix B. 
The procedvire is extremely tedious, and, in fact, vumecessary, since the 
blast load is not accurately knovn. Therefore, a more convenient method 
such as that established in Subsections 4.5.1 and 4.5*2 is desirable. 

In these subsections, the blast load acting on the arch has been 
idealized to tvo initially peaked triangular force pulses. For investigat- 
ing the effect of the blast load on the haunch of a I80 degree arch, the 
peak value of the first pulse is and the peak value of the second 

pulse is 0.06 p^. For similar investigation of a 120 degree arch the peak 
value of the first pulse is 0.45 vhereas the peak value of the second 
pulse is 0.04 P^. In each instance the duration of the first pulse is the 
time required for the pressure vave to engulf the arch, and the duration 
of the second pulse is equal to the effective duration of the bleist loading. 



TABLE V 

HARMONIC ANALYSIS OF BLAST LOADING ON l80° ARCH 
(p ('t) obteiined from Table B5> Appendix B, for t * 10 ms) 
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PEAK VALUE OF ANALOGOIB LOADS p (t) AMD p (t) FOUND BY HARMONIC ANALYSIS 
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Figure 4.1 Geometry of Blast Loading on a l80 Degree Central Angle 
Circular Segment Arch 



'dV 
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Figxare 4,2 Overpressure as a Function of Time 



200 
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Angle <f> , degrees 
(b) 120 Degree Arch 

Figure 4.3 Blast Loading on l80 and 120 Degree Central Angle 150 Inch Radi vis Arches at Time t = 10 ms 



Peak Overpressure, psi. Peak Overpressure, psi. 
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Maximum Value Pg(t), psi. 




Maximum Value Pg(t), psi. 
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Figure 4.4 Peak Values of Analogous Loads: Anti symmetrical p„(t) and 
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Symmetrical p (t) for a l8o Centreil Angle Arch 
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Figure U .5 Peak Values of Analogous Loads; Anti symmetrical and 

Symmetrical p (t) for a 120 Degree Central Angle Arch 
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CHAPTER V 

PASSIVE EARTH PRESSURE LQADIMG 

If there vere no soil surrounding the arch, it would deflect into 
the dotted configuration shown in Figia:e 5*1 under the unsynimetrical blast 
loading. However, the soil restricts this deflection by passive earth 
pressiure and the inertia of the soil mass. To evaluate the passive earth 
pressure, the segment of the arch on the leeward side from the springing line 
to the crown is assxuned to act as a retaining wall. The soil volvune and 
the relative plane of slip are shown in Figure 

Under actvial conditions, the body of soil located on the windward 
side of the arch woxild most certeinly be affected by the blast pressure and 
would have to move downward and laterally into the arch to load it. The 
energy required to accelerate this mass woxild lower the energy level in the 
blast pressure. In addition^ after the surface pressure wave passes the 
arch, there would be a vertical pressure of varying intensity acting downward 
on the soil body on the leeward side of the arch. This would increase the 
shear strength of the soil, the passive earth pressxare, and the inertia 
force. Evaluation of these later effects would be extremely difficxilt and 
is omitted in this study. This resxilts in a conservative approach. 

5.1 Unit Passive Earth Pressure 

i Rankine's formula for the total horizontal passive earth pressxrre 

force P^ acting on the vertical plane AV with a height h in Figirre 5*1 is 

P 

1 + sjLn 

where K = t rj — ^ , <P' being the angle of internal 



and 



w 



unit weight of the soil. 
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The unit horizontal pressxare eilong the plane AV in Figure 5*1 
varies directly as y varies, y "being the distance helov the grovtnd surface. 
This pressure then may be e^tpressed as 

Pp ® vyK (25) 

5.1.1 180 Degree Arches 

For a semicircular arch with no cover on the crown and <p' = 50 
degrees, since y = R(l - sin<p). Equation (25) becomes 

= wR(l - sincp) l qJ = 5vR(l - sin<p) (26) 

The angle <p is measiired from the horizontal radial line \diich joins the 
abutment with the center of the circular segment. 

In addition to the horizontal passive pressure, it is necessary 
to consider the weight of the soil above the arch rib. iMs soil exerts a 
unit pressure 

p^ = wR(l - sin<p) (27) 

Using Equations (l), (26), and (27), the resultant radial earth presaure is 

Pj.(^) = i (Py + Pp) + i (Pp - P^) cos 2<P 

= 2 ['vfi (1 - sin<p) + 3'wR (1 ” sin9) ] 

+ 5 [5'JR (1 “ sin<p) - wR (1 - sincp) ) cos 2<p 
= 2wR (1 - sin<p) + wR (1 - ain<p) cos 2<p (28) 

' In Figure 5.I, the angle p, which defines the inclination of the 

plane of shear AS with the veirtical plane AV, is 

p = ^ (90° + 9')> being the angle of internal friction 

of the soil 



i (90° + 50°) = 60° for <P’ = 30° 



!tt 
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From the geometry of Figure 5.1, the total veight of the soil 
associated vith the passive earth pressure is 

Wg = ^ R(R tan 6o) + R^ - •— J ” 1.08 til (29) 

Equation (28) is for an arch vith no cover on the crovn. If 
there is a depth of cover on the crovn, a uniformly distributed hori- 
zontal pressure 5'4I^/ as derived from Equation (25), and a uniformly dis- 
tributed vertical pressure must be added to the radial pressiire acting 

on the arch. For an arch with a centreil angle of l80 degrees, the resulting 
radial pressure due to H^, as obtained by Equation (l), is 

Pj.('P) = 5 (Py + Pp) + i (Pp - P^) cos 2<p 

= j (>fi + 5^ ) + 2 (5'ffl - vH ) cos 29 
0 O' 0 o 

* 2 vfl + vH cos 29 (50) 

c c ' ' 

The total radial pressxare on a l80 degree arch with a depth of cover 

is found by adding the pressures derived fl*om Equations (28) and (50). 

P^(*p) = 2 vR (1 - sin9) + vR (1 - sin9) cos 29 

+ 2vfl + vH cos 29 
c c 

5.1.2 120 Degree Arches 

Equation (28) is valid only for a semi-circular arch vith the 
angle 9 measured from the horizontal radial line. For an arch with no 
cover on the crovn and a central angle of 120 degrees. Equation (28) is 
modified by substituting (9 + 50*^) for 9. This yields 

p^(9) = 2\« (1 - sin (9 + 50 )] + vR (l - sin(9 + 50)] cos 2(9 + 50) 

= -wR ( 2 - ^ sin9 - ^ cos9 - ^ cos 59 

"II 1 

” 2 29 + ^ cos 29 ) 



(51) 
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If there is a depth of cover on the crown of the 120 degree arch, the 
resulting radial pressure due solely to is 

Pj.(9) = i (Pjj + P^) + i (Pjj - P^) cos 2cp 

= i (5vH + vH ) + i - wH ) cos 2(9 + 5°) 

c c c c 

= 2 wH + i ^ cos 29 - ^ sin 29 (52) 

C C d c 

• 

The angle 9 is meeisxired from the radial line Joining the abutment with the 
center . 

The total unit radiaJ. jjressure of a 120 degree arch having a 
depth of cover on the crown is found by adding Equations (5I) and (52). 
The resvilt is 

Pp(9) = wR ( 2 - sin9 - ^ cos9 - ^ cos 59 

- sin 29 + 5 cos 29) + 2 vffl 
d o 

+ 7 cos 29 - ^ wH sin 29 (35) 

C d c 

The total weight of the soil associated with this force is 

(R + J w 

= 0.542 wR^ + 1.T52 wRH^ + 0.866 wH^ 

c c 

5.2 Distributed Earth Pressure Load 

In considering the unsymmetrical blast loading on the arch rib, 
it is advantageous to resolve this loading into harmonic components since 
the resulting distributed loads ^ easier to treat 

analytically than the actual irregular blast load. This advantage does not 
exist for the passive earth pressiure loading, since original anadytical 
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e^^esslons are available which may be studied in detennining the effect 
of the passive earth pressure on the arch. However, resolution of the 
passive earth pressiare into harmonic loadings does feicilitate visualization 
of the faot that the earth pressure may be thought of as a distributed 

p^ sin in the arch rib. In addition, a desirable uni- 
formity of approach is provided for the entire problem. 

Equation (28) expresses the xinit radial pressure acting on a 
180 degree arch with zero depth of cover on the crown. The toteil force 
acting on the arch rib may be resolved into harmonic components by 
Equation ( 9 ) • 

= 2wR(1 - sin9) + wR(l - sin?) cos 2<P 

--hj ^ 

It 

p^ = [2wR(1 - sin?) + wR(1 - sin?) cos 2?]sin Rd? 

= 0.948 wR 
a 

P 2 = ^ = Q-509 «R 

% 

P5 “ ^ wR 

Equation ( 55 ) expresses the radial pressure acting on a 120 degree 
arch having a depth of cover on the crown. The total radial force re- 
svilting from this pressure may be resolved into harmonic loading components 



as follows: 
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p^(cp) = vR (2 - sin 9 - ^ coscp - i cob 59 

- ^ sin 29 + ^ cos 29) 

+ \(H (2 + I eos 2<P - ^ sin 2?) 

5n ' A / ^ 

« 

A 

[v/^ ^ ^ 



- ^ sin 29 + ^ cos 29] sin ^ Rd 9 



+ \ffl 



(2 + ^ cos 29 - ^ sin 29) sin ^9 Rd 9 j 



= O.m vR + 0 . 86 If \dl 



= ^ / ^ P (<P) sin 59 Rd 9 = O.I8I MR + 0.987 MR 

4«R '^o ^ ^ 



A 

= ~ / 5 p^( 9 ) sin §9 Rd 9 = O.I98 vR + 0-557 
|j«R '^o ^ ^ 



5 .5 Analogous Earth Presstire Load 

In order that the effect of the unsynnnetrical blast load on the 
arch rib may be determined conveniently, the blast load is represented as 
an analogous uniformly distributed radial load acting in con^aression on the 
windward side of the arch and in tension on the leeward side. This greatly 
facilitates the study of the dynamic effects of this load since the problem 
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of dealing with a time -dependent imsymmetrical load is reduced to that of 
dealing with a time -dependent uniform load. The problem may be simplified 
further by letting this uniform load define the penrameters of the dynamic 
load acting on the single-degree-of- freedom system. 

It is logical to extend the concept of analogous loads to the 
unsyrametrical passive earth pressure acting on the arch. In the first place, 
a desirable uniformity of approach to the entire problem is obtained. In 
addition, by defining the earth pressure load £is an anti symmetrical radial 
lead acting in ccxnpression on the leeward side of the arch and in tension 
on the windward side, the value of this pressure can be considered as inte- 
gral part of the resistance of the arch rib, since the direction of the 
earth pressure is directly opposite to that of the blast pressure. 
Furthermore, this resistance provided to the arch by the soil may be used 
to define the resistance of the single-degree-of- freedom system. 

The loading components ^ p^ sin defined in Section 5*2 , 
coupled with the moment coefficients from Table IV, are used to define the 
analogous loads for l80 and 120 degree arches. The resxilting analogous loads 
define the maximum, or yield, resistance given the arch rib by the soil, 
since the passive earth pressure is evaluated at its peak value. The 
resistance is evaluated eis a function of the moment in the arch rib at a 
point located by d/h, the same point used to define the analogous load for 
the blast pressure. 

For a l80 degree arch with no cover on the crown, the resisting 
moment due to the soil is 



M = -0.025(0.948 wR)R^ + 0 . 555 ( 0.509 vR)R^ + 0.088(l.l80 wR)R^ 
= 0.252 wR^ 



i]Q* 
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The analogous load Is 

0.252 wR^ 

P« = T* 

0.414 

= 0.611 vR (54) 

For a 120 degree arch with a depth of cover on the crown, the 
resisting moment due to the soil presstire is 

M = -o.oo9(o.iUwR + o.864'iai^)R^ + 0. 125(0. iSlwR + 0.987 wH^)R^ 

+ 0.050(0.198^® + 0 . 557 WH^)R^ 

= 0.054wR^ + 0.144^^R^ 

The analogous load is 

0.054wR^ + 0.144^ R^ 

P = 5 

0.155R^ 

=> 0.219WR + 0.929WH^ 



( 55 ) 
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Figure 5*1 Passive Eaz*th Pressure on a l8o Degree CeatreuL Angle Arch 
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CHAPTER VI 

ANALYSIS ^ ^ RESPONSE OF TM ARCH TO ^ DYNAMIC LOAD 

The arch rib is idealized as a single-degree-of-freedom system, 
and the flexural resistance and mass of the arch rib per se are ignored. A 
method of dynamic analysis developed by Dr. N. M. Nevmark, References (8) 
and (9)> is used to predict the response of the arch rib to the given 
pressure pulse. 

For the arches investigated in this thesis, it is found that the 
response of each is in the plastic, or yielded, part of the assumed elasto- 
plastic resistance diagram for the soil. Had the response been less than 
the deflection required to establish the peissive earth pressure, the dif- 
ferentieO. equation of motion stated in the introductory chapter of this study 
would have been used to find the maximum response. 



6.1 Method of Analysis 

For an initially peaked triangular force piilse acting on a single- 
degree-of-freedom system having an elasto-plastic resistance fvmction. 

Figure 6.1, Dr. Newmark presented in Reference (8) the following approximate 
equation to predict the maximum response of the structure; 



^max 




+ 




T 

1 + 0.7 ^ 
^1 



( 56 ) 



where 



^max ~ value of the pulse loading 

q^ = yield resistance of the structural element 

|i = ductility ratio = maximum deflection/yield deflection 




T 



natural period of vibration of the stiructiire = 2« 



= duration of the loading 






Another technique developed by Dr. Newmark, Reference ( 9 )^ can 
be used to define the dynamic behavior of the structure for other pulse 
shapes >dilch may be represented by several initially peaked triangular force 
pulses. Figure 6.*1. The relationship betveen the peak pressure and the 
marl mum response of the structure is defined approximately as 




vtoere F^, Fg = damage pressure levels for initially peaked 

triangular force pulses of duration t^/T and 
tg/T, respectively, acting separately 

c^, Cg = ratio of the peak pressure of the first and 

and second replacement triangles to the total 
peak pressure 

^1" ^2 “ respectively, the 

maximum value of the individual force pulses 



To apply the above described methods of dynamic analysis to the 
situation at hand it is assumed that: 

( 1 ) The configuration of the arch rib is defined by the radial 
deflection of the haunch at a point located at e/4. The structure is 
idealized as a single-degree-of-freedom system. 

( 2 ) The blast loading on the arch is considered to be a uniformly 
distributed anti symmetrical loading. The unit value of this loadli^g defines 
the forcing function acting on the single-degree-of-freedom system. 

(5) The shape of the blast loading force pulse consists of two 
initially peaked triangular force pulses. The duration of the first pulse, 
t^, is defined as the time required for the blast wave to engulf the arch 



72 



rib and rise to its maximum pressure at the leevard abutment^ a time equal 

to t plus t . The total dxiration of both pulses is defined as the effective 
a r 

duration of the blast loading, t^. 

(4) The resistance provided to the arch rib by the passive earth 
pressure is considered to be a uniformly distributed anti symmetrical load. 

The unit value of this resistance defines the resistance function of the 
single-degree-of-fheedom system. 

(5) The vel^t of the soil aissoclated with the passive earth 
pressvire is uniformly distributed along the arch rib. This vtnlt veight 
defines the vinit mass of the single->degree>of-freedom system. 

(6) The passive earth pressure reaches its maxi mum valize ^en 
the radial deflection of the arch at a point d/4 is 0.1^. The resistance 
diagram of this pressure is elasto-plastlc and defines the resistance func- 
tion of the idealized single-degree-of-freedom system. 

(7) The period of vibration of the arch is defined as 

T = 2ti \/mx^/^ , -vdiere m results from the unit veight found in (5)» is 
0.1^, and is the unit peussive earth pressure resistance found in (4). 

The general nature of assunqptlons (l) throvigh (7) above is shovn 
in Figure 6.2. 

Numerical values of the parameters in Equation (56) are presented 
in Table VIII. Using these parameters, the maximum responses of arches 
■with various radii subjected to ICX) and 200 psl overpressures have been 
de'bermined and are plot'bed in Figure 6.^. It is assvuoed that 'the arch is 
burled in soil with a seismic -velocity of 2000 fps* The resistance and 
mass of the arch rib have not been included, and soil arching effects are 



not considered 
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6.2 Sanple Calculation 

A sample calculation of the resxxsnse of an arch under an assumed 
set of conditions follows. 



Let: 



Cg = 2000 fps 

= 2 ft/ms 
w = 120* pcf 




W = 1 MI 
R * 200 in. 

» 16.67 ft. 

H = 0 
c 



P » 100 psi U = 5000 fps 



= 5 ft/ms 
6 - 180 degrees 



The parameters are: 

W * ^-57 = 0.57(100) = 57 psi 

p' = 0.06 P = 0.06(100) = 6 psi 
•*^max m \ / x' 

^ _ R - R 6in(<p - p) R - R sin9 

^1 “ C U 2C 

3 s 

16.67(1 - sin(0 - arctan |) ] 

= gJ 2 ^ + 

- 20 ms 

m 

= 400 ms 

T = 0.0076 n/r = 0.0076 'iioo = 0.107 sec. 

= 107 ms 

Xg = 0.001 R = 0.001(200) = 0.2 in. 

= 0.611 wR = 200 = 8.5 psi 



7 ^ 



The force pulse is shown in Figure 6.2. Equations (56) and (57) 
must be evaluated by trial and error. Reference (9)> Chart I, may also be 
used to evaluate Equation (56) . 

Assume that n = 12. For the first initially peaked triangular 

force pulse having a dxuration of 20 ms, the damage pressxure level necessary 

for a maximum deflection of x » 12x is 

m e 



P, 






max T fZ T 
% " 2|i - 1 + 1 + 0.7(T/t^) 

= ^ v/24 - 1 + 



^ ■ k 



20 * 



1 + 0.7(107/20) 



= 8.4 



For the second initially peaked triangular force pulse having a 
duration of 400 ms, the damage pressure level necessary for a ma-yitmim de- 
flection of X = 12x is 
m e 



°^max 

% 



10 



V24 - 1 + 



1 + 



0.7(107/400) 



1.2 



The small increaae in the peak value of the second pulse due to 

the extension of the pressure line from t=20mstot=0 ms, an increase 
20 

of (6) ss 0.5 psi, is ignored. 

The actual peak pressxure in the first force pulse is 31 psi; the 
actual peak pressure in the second force pulse is 6 psi; and the actual 
ratio of the total peak force to the resistance is 



^max _ 57 

% " 8.5 



4.57 



Equation (57) niay now be evalxxated to determine if the assumed 



value of |i = 12 is satisfactory. 



I 
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1.0 



■where 




0.84 




o.l6 



4.57 ( 2^ * f#) 



1.02 = - 1.0 



Therefore the aissuined value of 12 for n is satisfactory, and the 
maximum deflection is 

12 ( 0 . 2 ) 

2,4 Inches 



X = 

m 



6,5 Discxxssion of Resvilts 

Figvure 6.5 contains the quantitative results of the analyses of 
the response of various sized arches. The results may be qualitatively 
svunmarized as follows: 

(l) For any given overpressure level, the response of 

the arch increases as the raditis of the arch decreases. 

Biis is partially expledned by considering how the load 
on the arch and the resistance given to the arch by the 
earth pressure are defined. The bleist pressure is de- 
fined as a function of the peak overpressure only. The 
value of the axialogous, uniformly-distributed, anti- 
symmetrical load is not a function of the radius of the 
arch. The earth pressure is defined as a function of 
the radius of the arch and varies directly as the radius 
varies. Consequently, the response of the arch to any 



given overpressure varies indirectly as the radius 
of the arch varies. 

Other factors •vdiich probably influence the re- 
sponse of the arch have not been studied in this thesis. 

One of these is the value of the resistance of the arch 
per se. For arches of smaller radii> this resistance 
may be relatively large compared to the resistance given 
the arch by the earth pressxire. 

For equal radivis arches^ the response of the eurch 
to any given overpressure level increases as the central 
angle of the arch decreeises. 

The effect of the second triangular force pulse is 
very inportant since it has a relatively long duration. 
Several analyses of the response of a 120 degree eirch to 
a 100 psi overpressiare vere made, in ^diich the effect of 
this second force pulse vas not considered. The response 
of a 200 inch radius arch vas foxmd to be approximately 
three inches, con 5 >ared to a value of sixteen inches idien 
the effect of the second pulse is included. 

Increasing the depth of cover on the crown reduces 
the effect of the blast load on the arch. The effect of 
the depth of burial of the arch may be related conveniently 
to the ratio of the depth of burial of the arch abutment to 
one-half of the arch span. When this ratio is one, for 
arches having radii of 200 inches or more and central angles 
greater than 120 degrees, the flexural action caxised by an 
overpressxire of 100 psi resvilting from a 1 MP weapon may 



be ignored. For a 180 degree arch^ this ratio is one 
•vdien there is no cover on the crown, and the deflec- 
tion is approximately two inches. For a 120 degree 
arch, this ratio is approximately 0.9 "vdien there is 
a depth of cover of O.JR on the crown, and the deflec- 
tion is approximately one inch. 



I 





78 



TABLE VIII 

PARAMETERS FOR DYHAMIC AMALISIS OF SOIL RESISTAMCE 



Parameter 


l80 Degree Arch 

H = 0 
c 


120 Degree Arch 

H = 0 
c 


120 Degree Arch 

H = O. 5 R 
c 


P 

^max 




O.U5P^ 


0.45P^ 

m 


p’ ^ 


"0.06P 


0.04P 


0.04P 




m 


m 


m 


t. 


t + t 


t + t 


t + t 


1 


a r 


a r 


a r 




t. 




t. 


2 


1 


i 


i 


‘Je 


0.611vR 


0.219VR 


0 . 497 ^ 


X 


O.OOIR 


O.OOIR 


O.OOIR 


e 

m 


1.08wR^ 


0.342wR^ 


0.940\«^ 


« Rg 


I *Re 


1 * Rg 


T 


0.0076 n/r 


0.0087 ^/R 


0.0096 ^/R 



I 





79 



q(x) 



X 



Mass 



p(t) 



(a) Single-Degree-of -Freedom 
System 




(b) Initially Peaked 

Triangular Force Pulse 



q(x) 





(c) Elasto-Plastic Bilinear 
Resistance 



(d) Two Initially Peaked 
Triangular Force Pulses 



Figure 6.1 Idealized Single-Degree-of -Freedom System^ Force Pulses and 
Elasto-Plastic Bilinear Resistance 
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p(t) 





Figure 6.2 Idealized Loadiiig and Resistance Functions of an Underground 
Arch Rib Subjected to Blast Pressure 



I 




Deflection Of Haunch At -;p , inches 
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Figure 6.^ Eaunch Deflection vs Arch Radius 

W = 1 MT C„ » 2000 fps t = ^ . 

P_ as noted as noted 

01 c 



d as 33oted 
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CHAPTER VII 

QUALITATIVE DESIGN METHOD FOR ARCH RIBS 

The problem of the structural design of the arch rib has not been 
studied. Hovever, the following qualitative design procedxare, based on the 
methods used in this thesis, is suggested: 

(1) Idealize the surch as a single-degree-of -freedom system, the 
configxnration of -vdiich is defined by the radial deflection of the arch at 

e/4. 

(2) Define the fundamental period of vibration of the structure 
by the mass and resistance of the soil, and the deflection of the arch at 
e/4 necessajry to establish the full passive earth preesxire. 

(5) Define the forcing function ais consisting of two initially- 
peaked triangular force pulses, the first pulse having a dxuration of t + t , 
the second having a dvcration of t^. These forcing functions are imiformly- 
distributed, anti symmetrical, radial loads analogous to the actual blast 
loading. 

(4) Assume a ductility factor (i* The end objective is to 
design the arch rib, not to determine the response of the soil surrotmding 
the arch. However, since the arch cannot deflect without the soil ac- 
companying it, this ductility factor also defines the response of the soil. 
The assumed value of m should be larger than that desired for the actxial 
arch rib. 

(5) Evaluate the damage pressvire level equation to determine 
the total resistance required of the combined arch rib-soil system. 

(6) Subtract the soil resistance ftrom the total required 
resistance and design the arch rib with a resistance equal to the remainder. 
The usually recommended allowable stresses for dynamic design are to be used. 
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( 7 ) Determine the yield deflection at Q/k of the arch section 

ft 

designed in (6) . The yield deflection must be stated in the same munerical 
terms as those used in (2) for finding the period of vibration of the 
structvjre . 

(8) Combine the resistance of the soil and arch rib into an 
elasto-plastic resistance diagram. This is accomplished by making the area 
of the idealized bilinear resistance diagram equal to the area of the 
actual trilinear resistance diagram vhen both diagrams extend to a numerical- 
ly specified deflection. The ratio x is assxamed; the initial slope of 
the elastic part of the actual diagram is maintained; the value of is 
adjusted to make the areas equal. Determine a new period of vibration for 
the combined arch-rib soil system, including the resistance emd mass of the 
arch section designed in (6). 

( 9 ) Using the elasto-plastic resistance function found in (8), 
evaluate the damage pre8S\u*e level equation to determine the maximum 
response of the system. The resulting ratio of the maxi mum deflection to 
the yield deflection of the system is fovind and coirpared with the ratio 
x^Xg used in combining the soil and arch resistances. If these ratios 
are equal, the arch section designed in (6) is satisfeictory. If they are 
not equal, the arch section must be redesigned and (7) throxigh (9) repeated. 

The anticipated results of the design procedure outlined above 
are discussed below without proof. 

There are two principle ix>des of response of an arch to blast 
loading. Reference ( 10) . One of these is symmetrical^ corresponding to 

the compression in the arch, and the other is antisymmetrical, corresponding 
to the flexural deformation of the arch. Consider an underground arch and 
assume that the arch rib has been designed for the static load of the soil 
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cover and the compression load of the blast force. For the arch sizes, 
overpressures, emd weapon yield studied in this thesis, it is reasonable 
to expect that such arches woxild have a yield deflection on the haunch of 
at least one-half inch. For 100 psi overpressure, the tnaxlimun deflection 
e:^cted of the haunch, considering the resistance of the soil only, is 
two inches for l8o degree arches with zero depth of cover on the crown, and 
approximately one inch for 120 degree arches with O.JR cover on the crown. 
Therefore the ratio of the maximum deflection on the haunch to the allow- 
able elastic deflection of the arch rib does not exceed 4:1, even though 
the reduction of the maximum deflection due to energy absorbed by elastic 
or plastic deformation of the arch rib is disregarded. 

For moderately ductile structures, the above ratio, ductility 
factor u, ranges from 10 to ^ 0 , Reference (8). For designing flexural 
members, including reinforced concrete members, the ductility factor usually 
is assumed to be 20. 

It is therefore concluded that if the arch rib is designed so that 
it has a ductility factor of 20 for flexural action, and is covered as de- 
scribed above, no increase in the required resistance for static and blast 
compression load is necessary to accommodate the anti symmetrical mode of 
response for an overpressure of 100 psi. 
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CHAPTER VIII 

COMCLUSIOHS AHD RK^OMMEHDATIOMS 

8.1 Conclusions 

The fo3J.owing conclusions are drawn from the results obtained 
and the discussions thereof in the preceding chapters. 

(1) Harmonic analysis ixrovides a convenient means of analyzing 
unsymmetr i cal loads acting on arches. 

( 2 ) The peak values of dynamic loads euad their load- time 
curves may be determined by harmonically analyzing the load at suitable 
time intervals. 

( 3 ) The effects of the unsymmetr ical blast load acting on the 
arch may be stxidied conveniently through the use of analogous loads. When 
studying the effect of "Uie blast load on the havmch, the analogous load 
en 5 )loyed is a uniformly distributed radial load ■vdiich acts in con^iression 
on the windward side of the arch and in tension on the leeward side. The 
peak value of this load is approximately forty percent of the peak value 
of the overpressure. When studying the effect of the blast load on the 
crown, the analogous load employed is a uniformly distributed radial load 
•vdiich acts in compression on the central third of the arch and in tension 
on the outer thirds. The peak value of this load is approximately forty 
percent of the max1,mvim overpressure. 

The effect of the unsymmetriceil passive earth pressure load 
on the arch also may be studied throvtgh the use of an analogous load. This 
analogous load is a uniformly distributed radial load acting in compression 
on the leeward side of the arch and in tension on the windward side. It 
directly opposes the blast pressure load. Its value varies directly with 



the radius of the arch. 



■e 






I 
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(5) The peak value of the analogous load en 5 >loyed in studying 
the blast effects on the haunch of the arch Increases as the seismic 
velocity of the soil increases. Bie peak value of the euialogous load vised 
in stvdying the crovn of the arch decreases as the seismic velocity of the 
soil increases. 

(6) The relative severity of the dynamic load varies inversely 
with the ratio of the horizontal to vertical pressure in the soil. It is 
to be noted, however, that as this ratio increeises, there may be a tendency 
for the crown of the arch to move outward. 

(7) The relative rise time of the blast pressvire to its peak 
value has little effect on the peak value of the analogous loads. 

(8) The passive earth pressvire contributes materially to the re- 
sistance of the arch. An overpressvire of 100 povinds per square inch from a 
one megaton weapon causes a maximum response of approximately two inches 

in the l8o degree arches studied. Subjected to the same overpressvire, the 
response of a 120 degree arch with a depth of cover on the crown equal to 
three-tenths of the radius is approximately one inch. This response is 
considered negligible. 

( 9 ) The passive earth pressvire developed by the arch decreases as 
the central angle of the arch decreases. It is doubtful if arches with 
central singles of much less than 120 degrees would develop significant 
passive earth pressure under blast loading. 

( 10 ) The rigidity of an arch interferes with the deflection of 
the arch required to develop the passive earth pressure. This sviggests that 
flexible arches have certain advantages in blast resistant construction. 

( 11 ) The method of approach presented in this thesis may be 
extended to consider soils with different angles of internal friction. 
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different xmlt weights, and varying water table levels, as well m to 
cohesive soils, provided the cohesion of the soil can be stated in terms 
of an equivalent angle of internal friction. 

8.2 Recommendations 

It is recommended that the studies made in this thesis be 
extended to; 

ft 

(1) Consider how the blaist loading on the arch is affected by 
the soil on the windward side of the arch. 

(2) Consider how the passive earth pressture loading on the arch 
is affected by the Increased shear strength of the soil due to the veortical 
blast pressxire acting on the leeward side of the arch . 

(5) Coa^iare the arch behavior predicted in this thesis with 
available test data. 

{h) Develop a quantitative procedure for design of the arch rib 
under impulsive loading vdxich includes the effect of the passive earth 



• pressure 
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APPENDIX B 

NUMERICAL CALCULATION RESULTS 

This Appendix contains the results of the numerical harmonic 
analyses made. 

In Tables B1 through BIO are tabulated the vertical pressure-time 

values of blast pressure on a l80 degree arch vith zero depth of cover on 

the crown for various combinations of the variables stated in Chapter II. 

Since the 120 degree arch has the same radius, this table also provides the 

vertical pressure-time values for the blast loading on a 120 degree arch by 

using the values between JO and IJO degrees. 

Tables Bll and B12 contain the values of the coefficients p (t) 

n 

for the harmonic components of loading, as found by nxunerical analysis. 

These values were found by using the vertical pressures in Tables B1 through 
BIO, applying the coefficients from Table I to convert the vertical pressure 
to radial pressinre, and then harmonically analyzing these radial pressures. 
These tables also contain the values of the ansdogous loads p„(t) and p (t). 
The following sign convention is applicable to these loads: 

(s-) Pn(^) positive vden the load suits inwardly 

(b) p (t) is positive vhen the uniformly distributed load 
a§ts in compression on the windward half of the arch 
and in tension on the leeward side 

(c) p (t) is positive when the uniformly distributed load 
aSts in compression on the central third of the arch 
and in tension on the outer thirds. 

Originally it was intended to study only twenty- five foot radius 

arches, some with a rise time equal to one-half the transit time and the 

remainder with a rise time equed to the transit time. However, after nvimerous 

calculations had been made, an error amovmting to a multiple of two was found. 
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In order to vise the previous vork it became necessary to decrease the arch 
radius to 12.5 feet, and to increase the rise times of some of the arches 
to equal the transit time, and of the remainder to twice the transit time. 
This accounts for the data shoving t^ = 2t^. After the error was dis- 
covered, only arches with a twenty-five foot radius and t = were 

r Xf 

studied. The results are of such a nature that general conclusions may be 
drawn without regaJrd to the disparity in the arch sizes. 

Figures B1 through B6 are typical graphs of the values of the 
coefficient p (t) and the analogous loads p (t) 6ind p (t). 
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Table BJ. VerticeO. pressure-time plot for l80 degree 150 inch radius arch subjected to 
100 psi overpressure from 1 MT veajwn with C = 2000 fps and t = t . 
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Table B4. Vertical preBsure-time plot for l80 degree 150 inch radius arch subjected to 
100 psi overpressure from 1 ME weapon with C = 2000 fps and t = 2t. 
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Table B7- Verticeil pressure-time plot for l8o degree 15O inch radius arch subjected to 
200 psi overj)res8\ire from 1 MT weapon with C = 2000 fps and t = t. 
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Table B8. Vertical pressure-time plot for l8o degree 150 inch radius arch subjected to 
200 psi overpressiire from 1 MP weapon with C = 2000 fps euad t = 2t. 
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Table fill. Values of the coefficients p (t) for the harmonic components of loading p^(t) sin , 

n n C7 

anti symmetrical p (t) and symmetrical p (t) analogous uniform loadings, and the uniform 
a s 

all-around con^jresslon p (t). (All loadings are in pounds per square inch. Time is in 
milliseconds. R = I 50 inches. W = 1 MT. Other variables are as noted.) 
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Table B12. Values of the coefficients p (t) for the harmonic components of leading p (t) sin — ^ , 

n n c/ 

antlsyminetrlcal p (t) and symmetrlcaJ. p (t) euoalogous uniform loadings^ and the uniform 
a s 

all-around con^xression p^(t). (All loadings are in pounds per square inch. Time is in 
milliseconds. R = 5 OO inches. W * 1 ffT. Other variables are as noted.) 
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Figure b 4- Antisynnnetrical p (t) and Symmetrical p (t) Analogous Loads 

SL S 

with R = 150 inches, = 2000 fps, = 200 psi, t^ = t^, 

0 = 180 degrees, W = 1 ME, and noted 
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